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1 Introduction 

Unaveraged multichannel MEG measurements ex¬ 
hibit spatial and temporal correlations that can be 
used for quantitative data interpretation. One method 
for such interpretation is based on linearly con¬ 
strained minimum variance beamformers and pro¬ 
vides information about the 3D distribution of the 
source power. We shall consider a linear beamfor- 
mer (LB) and a nonlinear beamformer, called Syn¬ 
thetic Aperture Magnetometry (SAM). We shall 
show that under simplifying assumptions, it is pos¬ 
sible to describe the beamformer outputs analytical¬ 
ly. The analytic results will be compared with simu¬ 
lations and it will be demonstrated that in most ca¬ 
ses, the SAM produces sharper images than the LB. 
The dipole moment at a specific location 0 in the 
brain can be decomposed into orthogonal compo¬ 
nents. We shall assume a simplified case of MEG 
sources in a spherical conductor where only two 
tangential components are relevant (the extension to 
three components is straightforward). Two signal 
space vectors B ek , k = 1, 2 corresponding to the two 
dipole components can be combined into a gain ma¬ 
trix H e = [B 0 i, B e2 ]. Let the vector m(t) represent the 
instantaneous MEG measurement. Define a spatial 
filter centered on a location 0 as y(t) = W^m(t), 
where W 9 is the weight matrix. The elements of m(t) 
are assumed to be zero mean stationary processes 
and the beamformer weights are determined by 
minimizing the projected power Pe = WqCW 0 
subject to the condition W^Hq =1, where C is the 
covariance matrix of the measurement and I is the 
identity matrix. The weights are given by [1]: 

W 0 =(hTc- 1 H 9 )' 1 HTc- 1 (1) 

The weights can be used to compute various pa¬ 
rameters corresponding to the location 0, e.g., the 
time course of the dipole magnitude variation or the 
local power, or an especially useful quantity, the 
power normalized by the noise projected by the 
beamformer [1], This quantity will be called a 
pseudo-Z in this paper [2], A spatial image of the 
brain activity can be obtained by dividing the brain 
into voxels and directing the beamformer to the in¬ 


dividual voxels one at a time (a 3D search). 

The procedure where the dipole moment at each lo¬ 
cation is decomposed into two (or three) orthogonal 
components corresponds to the LB. An alternate ap¬ 
proach is to construct a nonlinear beamformer, 
where at each 3D location, the beamformer is syn¬ 
thesized only for one current dipole direction. This 
direction is rotated in the tangential plane and the 
orientation at which some beamformer parameter is 
maximized (e.g. the projected power or the pseudo- 
Z) corresponds to a solution for that location. This 
search is 4-dimensional and the procedure corre¬ 
sponds to the SAM. By retaining only one forward 
solution vector, H 0 = B 0 , the SAM weights are ob¬ 
tained from Eq.l as We = b£C _1 /(bJC _1 B 0 ). The 
location indicator 0 now represents three Cartesian 
coordinates plus the angle in the tangential plane. 

2 Methods 

First, the analytic derivation of the signal space 
vector parameters (norms and angles) will be out¬ 
lined, followed by the derivation of the analytic ex¬ 
pressions for the SAM and the LB. 

The geometry for the analytic calculations is shown 
in Fig.l. The sensor coverage extends from the ver¬ 
tex to a maximum angle (3 max and the sensor density 
is infinitely high. For simplicity, only radial mag¬ 
netometers are considered; however the treatment 
can be extended to radial gradiometers. A dipole q 0 
positioned at a distance a 0 from the model sphere 
center and oriented at an angle t|/ 0 relative to the xi 
axis will produce a field B 0 (r, <j), (3) at an observation 
point (r, <|), (3), where (|) is the azimuth and (3 is the 
declination. In spherical coordinates, this field is 
given by B 0 = b 0 f 0 sin(<|)-\|/ 0 ), where b 0 = |0. o q o /(4Jt) 
and f 0 = f(a 0 ,r,(3) = a 0 sinp/(ao+ r 2 -2a 0 rcos(3) 3/2 
[3], Similarly, the field due to a dipole at a probe lo¬ 
cation 0, oriented along t|/ 0 = \|/ 0 + A\|/ will be given 
by B 0 . The vector norms and the angle between the 
vectors B 0 and B 0 are computed in terms of their dot 
products as an integral over the sensor surface, 
B 0 B e = b2MI oe cosA\|//[2(l-cosP max )], where 
M is the number of channels and the integral Io 0 = 
Jp=o fofe sin(3d(3 depends only on a,,, a 0 , and r. The 




Figure 1: Geometry for the beamformer analysis. 
Sensor shell with radius r is concentric with the 
model sphere and is rotationally symmetric around 
the x 3 axis. A tangential dipole q a is positioned on 
the x 3 axis. The SAM output is probed along the x 3 
axis, the probe distance from the model sphere 
center is and its moment vector q e is oriented at 
an angle Ay/ relative to q 0 . The LB output is also 
probed along the x 3 axis (for clarity shown below the 
dipole q a ), and the probe has two orthogonal 
directions q e i and q g2 ■ The signal space vectors per 
unit dipole moment corresponding to the above 
dipoles and probes are B 0 , B g B gl , and B g2 . 

norms ||B 0 || and ||B e || are obtained by setting ‘O’ = 
‘0’ (which yield integrals loo and I 0e ). Angle X be¬ 
tween the vectors B 0 and B 0 is given by 

cos A = = I° 9 cos A\|/ = Ecos A\|/ (2) 

IIBoPell 

Eq.2 shows that in the axially symmetrical case, the 
angle between the two signal space vectors separates 
into the contribution due to the dipole orientation 
difference, Av|/, and the contribution E due to the di¬ 
pole and probe positions and their depth difference. 
If the probe was on the dipole, 0 -» O, a 0 = a„, then 
A = ±Ai|/, i.e. the angle between the signal space 
vectors is equal to the difference of their corre¬ 
sponding dipole orientations. If the dipole and probe 
were aligned, Av|/ = 0, then cosA = E. The depend¬ 
ence of E on the dipole positions and separation Ax 
resembles a cosine function and in the vicinity of the 
position O it can be expanded as E = 1 - yAx 2 , 
where y is a constant. 

The integrals Io 0 , loo, and I 00 (and therefore the E 
and |]B||) can be obtained analytically and they have 
an especially simple form for hemispherical (|3 max = 
7t/2) and spherical ((3 max = 7t) sensor shells. Complete 
analytical expressions for E, y, and ||B|| are not 
shown here. 


To derive analytic expressions for the beamformers, 
assume that there is only one dipole source, q 0 , em¬ 
bedded in a spherical conducting medium at position 
O, Fig. 1, and only the random sensor noise is acting 
on the system. For infinitely long integration, the 
noise covariance matrix is diagonal and if all sensors 
had identical noise, v, it would be given as E = v 2 I. 
The covariance matrix corresponding to one dipole 
q 0 and the random noise V is then C = V 2 I + q 2 BoBQ. 
This matrix can be inverted [4] to yield C' 1 =v' 2 (I - 
FB 0 Bq /|Bo)]\ where F = 1/[1 + v 2 /(q 2 ||B 0 || 2 )] < 1. 
The C' 1 and the weights shown in the introduction 
are used to obtain the beamformer parameters. 

First, consider the SAM and assume that the location 
0 is being probed while the dipole is at O, Fig.l. The 
projected power and the noise power are 
p e = W e T CW e and N e = W e T EW e . For simplicity, 
only the Z 2 = P 0 /N 0 will be discussed. Substituting 
the weights W 9 and the inverted covariance matrix 
C" 1 , the Z is a function of F and angle A between the 
signal space vectors B 0 and B 0 and is given by: 

2 B^C 'Be _ 1 - Fcos 2 A 

SAM v 2 BgC _1 C _1 B e 1-F(2-F)cos 2 A 

If no dipole was present, then q 0 = 0, F = 0 and the 
background value of Z S am would be obtained as 
Zsam bckgr = 1 . If the probe was on the dipole and ori¬ 
ented parallel to it, i.e. 0 —> O and A = 0, then 
Z 2 amo= 1 + qol|Bo|| 2 /v 2 . For small noise, 

v«q 0 ||B 0 || (large SNR), and the probe aligned 
with the dipole, Av|/ = 0, the half-power width of the 
Z|am depth profile peak can be determined using 
the expansion for E in the dipole vicinity. The first 
expansion term yields the peak width as wsam = 
v'2 / yV 2 /(qo||B 0 II 2 ). For large noise (small SNR), 
the SAM peak width approaches V 2(2-^2)/y . 

The analytic expression for the LB can be derived 
similarly. Since in this case the weights are matrices 
and not vectors, the projected power and the noise 
power are given as traces of the corresponding ma¬ 
trices, e.g. Pe= Tr{ (HQC _1 He)' 1 } [1]. After substi¬ 
tution of the inverted covariance matrix and labori¬ 
ous algebra, the pseudo-Z is obtained analytically. 
The general result is rather complicated, however, a 
significant simplification is realized by assuming 
axial symmetry. In this case Z 2 B simplifies to 


z2 _ 2-FS 2 (3-FS 2 ) 

LB 2- FS 2 (4 - F- FS 2 ) 


( 4 ) 















The Z LB depends only on the probe position through 
S and is independent of the probe orientation (be¬ 
cause two orthogonal tangential projections were 
used). The background value is again obtained for q 0 
= 0 as Z L b bckgr = 1. When the probe is on the dipole, 
then S = 1 and Z 2 LBo = 1 + q§||B 0 ]| 2 /(2v 2 ). If the 
noise is small (large SNR), then ZlBo~ Z 2 SAM c/25 
because the LB projects twice as large noise power 
as the SAM. The width of the LB depth peak for 
large SNR is w LB = (2 / Vy )V 2 /(q§ ||B 0 11 2 ) = V2 w SA m. 
For small SNR, the LB peak width is the same as 
that for the SAM. 

3 Results 

The accuracy of the analytic equations for the angle 
between the forward solutions is tested against 
simulations in Fig.2. In the immediate vicinity of the 
dipole, a e ~ a 0 and the analytic equations, the ap¬ 
proximate expansions, and simulations for the sen¬ 
sor arrays with interchannel spacing in the range 
from 1.3 to 6.5 cm agree well. Farther from the di¬ 
pole, the expansion results deviate from the simula¬ 
tions. The exact equations agree well with the 
simulations over a wide range of the probe positions 
when the channel separation is < 3 cm, indicating 
that the assumption of the infinite sensor density is 
not too severe. 

The expansion parameter y is shown in Fig.3.a as a 
function of the dipole distance from the model 
sphere center. For shallower dipoles, a 0 > 5 cm, y 
depends only weakly on the sensor coverage 
(hemispherical or spherical). For deep dipoles, the 
dependence of y on the coverage becomes strong. 
The ratio of the LB to the SAM peak widths is 



Figure 2: Angle between two forward solutions 
determined by analytic equations, expansion, and by 
simulations. The probe and dipole have identical 
orientations, Ay = 0. The simulations use discrete 
sensor arrays with different numbers of channels. 
Hemispherical sensor shell, radius r = 10.7 cm. 




Figure 3: Expansion parameter yand ratio of the LB 
and the SAM peak widths. Hemispherical and 
spherical sensor shells, r = 10.7 cm. (a) Magnitude 
of the coefficient y in the expansion E = 1 - yAx 2 ; 
(b) Analytically computed ratio of the LB to the 
SAM peak widths vs dipole magnitude, V = 3 JT 
rmsNHz, BW = 15 Hz, a 0 = 3, 5, and 7 cm. For a 0 = 
3 cm and q a < 0.7 nA.m the peak width cannot be 
determined because the background becomes larger 
than 1/2 peak magnitude. 

shown in Fig.3.b. For larger SNRs, the ratio 
approaches the expected V2 (for the dipole 
magnitudes of q 0 > 2 to 5 nA.m). For small SNR the 
ratio approaches 1. When the noise magnitude is 
changed, the curves shift to different qs. 

The analytic and simulated depth and angle profiles 
are compared in Fig.4 and in all cases the agreement 
is good. The depth profiles in Fig.4.a show that the 
LB peak magnitudes are 1/2 of that obtained with 
the SAM (note that a/z 2 -1 is plotted). Close to the 
peaks, the analytic results are slightly larger than the 
simulation results, the difference being smaller for 
larger channel densities. Also, increasing the chan¬ 
nel density reduces the peak width and increases the 
peak amplitude. The Z 2 as a function of angle in 
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Figure 4: Comparison of analytic results and 
simulations for the SAM and the LB. Hemispherical 
sensor shell, r = 10.7 cm, V = 3 JT rms/'JHz, BW = 
15 Hz, q a = 1 nA.m, a B = 5 cm. Simulations perfor¬ 
med with 100 and 316 channel sensor arrays (plots 
have different scales for different numbers of chan¬ 
nels). (a) Depth profiles when the probe is aligned 
with the dipole, Ay = 0; (b) Dependence on angle, 
a rotating probe is at the dipole location, a e = a„. 















Figure 5: Comparison of the analytically computed 
SAM peak widths with the approximate equations 
(expansions of E) and simulations. Hemispherical 
sensor shell, r = 10.7 cm, V = 3 fT rmsNHz, BW = 
15 Hz, q 0 = 1 nA.m, a 0 = 3, 5, and 7 cm, number of 
channels = 100. 

Fig.4.b shows that the LB output is independent of 
the probe orientation, while the SAM exhibits sharp 
peaks when the probe is aligned with the dipole. 
Similar to the depth profiles, the SAM peaks as a 
function of the angle are also sharper when the 
channel density is higher. 

Comparison of analytic and simulated peak widths 
for the SAM is shown in Fig.5. The analytic equa¬ 
tions predict the peak widths with excellent accu¬ 
racy. If the peak widths in Fig.5 were replotted as a 
function of Z 2 , they would become close to each 
other and could be described by a universal line. 
Simulated contours corresponding to the SAM and 
LB peaks are shown in Fig.6. The contours are 
elongated along the depth and as predicted 
analytically, the LB contour corresponding to 1/2 
Zp ea k is roughly V2 times broader than the SAM 
contour (in any direction). Figure 6 was computed 
for radial gradiometers and indicates that even 
though the analytic equations are strictly valid only 
for magnetometers, they also correctly predict radial 
gradiometer behaviour. 

4 Discussion 

MEG measurements are usually processed numeri¬ 
cally because the large number of channels makes 
the results analytically intractable. However, for the 
axially symmetric case, simple analytic equations 
for the SAM and the LB exist and can be used to 
predict their behavior without resorting to simula¬ 
tions. 

The LBs project twice as much noise as the SAM, 
their pseudo-Z 2 peak magnitudes are 1/2 of that ob¬ 
tained for the SAM, and their peaks are V2 wider 
than the SAM peaks. The peak widths of the two 
beamformers approach the same value for low SNR. 
Thus even though the SAM is computationally more 



Figure 6: Constant Z~ SAM and LB contours corre¬ 
sponding to a dipole in random noise. The Z 2 was 
normalized to 1 at the peak before the contours were 
constructed, dipole positioned at (xi, x 2 , xf) = (0, -2, 
5.66) cm, q a = 3 nA.m, yr = n/2, v = 3 JT/fHz, BW = 
100 Hz, 100 radial gradiometers with baseline b = 5 
cm, sensor shell extends from vertex to 1.66 rad, 
peak Z 2 sam ~ 14.3, peak Z 2 LB ~ 6.8. Solid contours 
correspond to half-peak values. The dotted line 
passes through the dipole position and the model 
sphere center, (a) SAM; (b) LB. 


intensive, its images are sharper. The present analy¬ 
sis considered only tangential dipoles with two 
components. If three dipole components were con¬ 
sidered, the noise projected by the LBs would re¬ 
duce their Z 2 peaks to 1/3 of the SAM peaks. 

The Z 2 peak width is inversely proportional to the 
number of channels, M. This monotonic decrease is 
a consequence of the random noise assumption (and 
is similar to the dipole fits which for random noise 
exhibit a scatter proportional to 1/VM). In the pres¬ 
ence of correlated noise, the reduction of the 
beamformer peak width with M reaches a saturation 
that depends on the random-to-correlated noise ratio 
and on the noise correlation distance. 
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